Entanglement between nitrogen vacancy spins in diamond controlled by a nano 

mechanical resonator 
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We suggest a new type of nano-electromechanical resonator, the functionality of which is based 
on a magnetic field induced deflection of an appropriate cantilever that oscillates between nitrogen 
vacancy (NV) spins in daimond. Specifically, we consider a Si(lOO) cantilever coated with a thin 
magnetic Ni film. Magnetoelastic stress and magnetic-field induced torque are utilized to induce a 
controlled cantilever deflection. It is shown that, depending on the value of the system parameters, 
the induced asymmetry of the cantilever deflection substantially modifles the characteristics of the 
system. In particular, the coupling strength between the NV spins and the degree of entanglement 
can be controlled through magnetoelastic stress and magnetic-fleld induced torque effects. Our 
theoretical proposal can be implemented experimentally with the potential of increasing several 
times the coupling strength between the NV spins as compared to the maximal coupling strength 
reported before in P. Rabl, et al. Phys. Rev. B 79, 041302(R) (2009). 



I. INTRODUCTION 

Nano-electromcchanical resonators (NEMs) arc at- 
tracting intense research efforts due to a number of fa- 
vorable properties such as the high sensitivity and the 
swift response to an external force with a low power con- 
sumption. This makes NEMs attractive for applications, 
e.g. for microwave switches, nano-mechanical memory 
elements, and for single molecule sensingir— . The role of 
quantized mechanical motion coupled to other quantum 
degrees of freedom as well as the influence of dissipation 
and noise are important issues in NEMs research with nu- 
merous findings and demonstrations of applications, e.g. 
Rcfs.— and further references therein. These studies 
also evidence the potential of NEMs for studying funda- 
mental questions concerning the quantum-classical inter- 
relation and issues related to entanglement and quantum 
correlations. Particularly interesting for the present work 
are spin states in nitrogen vacancy (NV) impurities in 
diamond^ir— as utilized for quantum information stud- 
ies. Their quantum mechanical properties can be mapped 
onto effective two-level systems which possess very long 
decoherence times even at a room temperature. On the 
other hand, they allow for a high degree of tunability via 
external magnetic fields which renders possible the use of 
NV impurity spins as sensors. For instance, a detection 
of a single electronic spin by a classical cantilever was 
demonstrated in2i. Magnetic tips attached on the free 
end of a cantilever generate a magnetic field gradient 
during oscillations, which induces a magnetic coupling 
between the nanomechanical oscillator and the spin sys- 
tem. Using the backaction due to this coupling one can 
readout the cantilever motio n^^'^^ . 

One further fascinating application for the NV spin 
based nano-electromcchanical resonator system is mag- 



netic resonance force microscopy (MRFM) which was 
proposed to improve the detection resolution for the 
three-dimensional imaging of macromolecules^. Oper- 
ation of the MRFM is based on the detection of the 
magnetic force between a ferromagnetic tip and spins. 
The possibility to achieve strong (up to 0.1 MHz or even 
stronger) coupling between the quantized motion of the 
nano-resonator and NV impurity spin was demonstrated 
in Refi^. However, an efficient control of the coupling 
strength and the development of ultrasensitive cantilever- 
based force sensors is still a fundamental challenge. Yet 
another application is based on using such a resonator 
in order to produce controllable entanglement between 
spins. Despite rather long decoherence time scales it has 
not been demonstrated yet, to our knowledge. 

In this paper we propose a new type of a nano- 
electromechanical resonator system, the functionality of 
which is based on the NV spin controlled by a magnetic 
field. The nano-resonator is a cantilever, which is covered 
with a magnetic film, such that we can exploit magnetoe- 
lastic stress or magnetic torque effects for a controllable 
cantilever defiection. At the free end of the cantilever 
magnetic tips are mounted, as indicated in Fig. 1. The 
magnetic field allows for a full control of the cantilever 
configuration via the amplitude of the applied magnetic 
field. This is the major novelty of the current proposal. 
As will be demonstrated below, the field-induced deflec- 
tion of the cantilever has a major impact on the inter- 
action strength between the NV spins and the magnetic 
tips, as well as on the strength of the indirect interaction 
between the spins mediated by the cantilever. We will 
show that the asymmetry of the nanomechanical system 
with respect to the shape of the cantilever controlled by 
magnetic fields drastically changes the degree of entan- 
glement. 
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The paper is organized as follows. In the next section 
we introduce our system as well as the model Hamilto- 
nian and discuss the relevant experimental parameters. 
Section lllll describes the magnetic field induced cantilever 
deflection. The subsequent Section IIVI focuses on the ef- 
fects of the induced indirect interaction of the impurities 
mediated by the cantilever coupling. The influence of 
the asymmetry coupling on the degree of entanglement is 
considered in Section |Vl Finally, Section I VII summarizes 
the results gives some perspectives for further progress. 



a) Spin 1 




FIG. 1. Schematics of the nanomechanical resonator system 
discussed in the present project, a) The distances /ii,2 be- 
tween the NV spins and magnetic tips can be adjusted by 
the field-induced defiection of the cantilever. The deflection 
is controlled by the applied magnetic field Bo. In the absence 
of the external constant magnetic field Bo = the system 
becomes symmetric, and /ii = /12 = fto ~ 25 nm b). Mi- 
crowaves are used to drive spin transitions at the NV centers 
of spin 1 and spin 2. b) Two magnetic tips are mounted on 
the free end of the cantilever, which mechanically oscillates 
in the z-axis direction. The length of the cantilever is of the 
order L ~ 3000 nm . The cantilever surface is covered with a 
10 nm thin Ni film. See Section Hill for further details. 



II. THEORETICAL MODEL 

The system of our interest is schematically shown in 
Fig. [TJ Magnetic tips are attached to the free ends of 
the cantilever on both sides. The deflection of the can- 
tilever and consequently the distance between magnetic 
tips and spins /ii_2 is controlled via the magnetic field Bq 
applied along the Z axis. The nitrogen vacancy center in 
diamond consists of a substitutional nitrogen atom with 




FIG. 2. Configuration of the energy levels. Q.± stand for the 
Rabi frequencies between the ground and the excited levels 
(I >,| -1 >) and (I >,| 1 >). 5± denote the detuning 
between microwave frequency cjo and transition frequencies. 

an adjacent vacancy. The total spin of the many-electron 
orbital ground state of the NV center is described by the 
spin triplet 5' = 1, tos = ~ljOil- States with differ- 
ent |ms| are separated by a zero field splitting barrier— 
which is of the order wq = 2.88 GHz. This splitting is an 
intrinsic property of the NV spin system and originates^ 
from the effect of the spin-orbit and spin-spin interactions 
leading to the single-axis spin anisotropy DS^ ~ Huiq. In 
the theoretical description we set ft = 1. The role of the 
applied external magnetic field is twofold: First of all, 
the applied external magnetic field (IbBq < fiwo, due to 
the Zeeman shift proportional to BqSz, removes the de- 
generacy of the levels | — 1 >, |1 >. Besides, as will be 
shown below, the external magnetic field, due to the thin 
magnetic Ni film deposited on the cantilever, modifies 
the shape of the cantilever. We will demonstrate that 
the asymmetry of the cantilever position between spins I 
and 2 (Fig. 1), induced by an external magnetic field, has 
important consequences for the coupling strength and the 
entanglement between the NV spins. Therefore, one can 
control the degree of entanglement and increase the inter- 
action strength between the NV spins. This constitutes 
an important advance in nanomechanics. 

The Hamiltonian of the single NV spin system reads^i 

^Nv = J2 {-^^ I*) + Y + 1*) • (1) 

In the case of a weak magnetic field i^bBq ^ fttJo 
(<C 30 mT) one can neglect level splitting and set (5_ = 
(5+, ri_ = fi+. In this case the Hamiltonian ([T]) cou- 
ples the ground state |0 > to the "bright" superposition 

of the excited states |6 >= — f(I ^ 1 > +|1 > )> while 

1 

the "dark" state \d >= —^{\ — 1>— |l>)is readout 

V2 

and decoupled from the process. Since only two states 
are involved, the NV spin triplet can be described via 
a S* = 1/2 pseudo-spin model. If the external constant 
magnetic field is strong enough, the splitting between 
the levels | — 1 >, |1 > is larger and the system becomes 
identical to the three level generalized Jaynes-Cummings 
model in the so called lambda configuration^2iS (See Fig. 
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A peculiarity of the generalized Jayncs-Cummings model 
is the expectation of two different transition frequencies 
between states (|0 >, | - 1 >) and (|0 >,|1 >). The 
Hamiltonian of the system holds a SU (3) symmetry and 
can be specified in terms of the Gell-Mann generators^. 
The three-level Jaynes-Cummings model is exactly solv- 
able in the general case. If the rf-field contains only 
one frequency resonant to the transition (|0 >, | — 1 >) 
the system can be reduced to the 5 = 1/2 pseudo-spin 
model which will describe transitions between the states 
I — 1 >, |0 > only, since the transition between the lev- 
els |1 >, |0 > is off resonance and therefore forbidden. 
However, the Rabi frequency f2_ of the transition be- 
tween the states | — 1 >, |0 > in this case is different 
from the transition frequency between the ground state 
|0 > and the bright superposition of the excited states 

\b >= ^^(1 — 1 > +|1 > )■ Nevertheless, in both cases 
v2 

our system can be described via the effective two-level 
model with a different Rabi transition frequency. The 
Hamiltonian of the system in the frame rotating with 
the frequency of the rf-field has the form 

Hs = 'S 1-1) (-11 + (o| + |0) (-11) . (2) 

Here S ~ ujq — fl{B) is the detuning between the mi- 
crowave frequency and the intrinsic frequency of the 
spins. The Rabi frequency of the transition between up 
and down spin states D,(B) depends on the amplitude of 
the magnetic field and, thus, includes both limiting cases. 
In the case of zero external constant magnetic field B = 
and zero splitting between the levels | — 1 >, |1 > the Rabi 
frequency for the transition between the bright \b > and 
the ground |0 > states is equal to ri(0) = Jlo, while in 
the case of a nonzero external field and nonzero splitting 
for the transition between levels | — 1 >, |0 > the Rabi 
frequency is equal to: 

n (B) = Oo - AO (B) , An (B) = fiB {Bo + . (3) 

Here Bq is the amplitude of the external constant mag- 
netic field applied on the system, i?i„s is the magnetic 
field produced by the magnetic film on the cantilever at 
the position of the spins Spin 1 and 2, and it is Bq > i?i„s- 
The eigenbasis of the Hamiltonian ([2]) is given by the fol- 
lowing states^^ 

|5) = cos (0/2) 1-1)+ sin (0/2) |0), 

|e) = - sin (9/2) |-1) + cos (6/2) |0) , (4) 

tanO = — r. 



In the basis (|3]) the components of the pseudo-spin oper- 
ator have the form 

= |e) (e| - 1.9) (5| , ^+ = \e) {g\ , a_ = \g) (e| (5) 

while the Z component of the spin Sz = 2(1'^ 

1 >< — 1|) and the Hamiltonian of the system Hs (given 



by eq. ([2])) reads 

5*2 = ^ [cos0(Tz + sin (tT+ + o-_)] , 

Hs = ^Lo<7,, cj = {n^ + 5')'^\ (6) 

Taking into account eqs. we have for the Hamil- 

tonian of a single NV spin interacting with the cantilever 

H = ^uja'' + a;,.(a+a + 1/2) + ^ (a+ + a) x 

[cosher"' + sin 6'(ct+ +0-")] . (7) 

Now we generalize the model given by eq. ([7]) for the 
system of two spins (see Fig. [1]). In the absence of an 
applied magnetic field the system is symmetric and the 
distance between the spins and the cantilever is equal 
hi = h2 = h. However, due to the magnetic field in- 
duced cantilever defiection, the external magnetic field 
leads to an asymmetry /ii 7^ /12. The imposed asymme- 
try Ah{Bo) = 2\hi - h2\l{hi + /12), < Ah{Bo) < 1 
can be quantified in terms of the amplitude of applied 
magnetic field Bq and the geometric and material char- 
acteristics of the cantilever, as outline din Section Hill 

The amplitude of zero point oscillations of the can- 
tilever is of the order of oq = ^/h/(2jnuj^ sa 5 x 10~-^'^?n, 
where m « 6x 10""'^'' kg is the resonator mass and 
uJi-/{2tt) w 4 MHz is the first resonance frequency of the 
cantilever. Details on the cantilever are given in Sec- 
tion 3. 

This oscillation amplitude is definitely smaller than the 
asymmetry imposed by the external magnetic field, i.e. 
Ah{Bo) 3> flQ. Nevertheless, oscillations of the cantilever 
produce a varying magnetic field due to the attached 
magnetic tips, which is proportional to the oscillation 
amplitude. The key consequence of the asymmetry is 
that constants of the interaction between the spins and 
the magnetic tips Ai^2 are different for different spins and 
the Rabi transition frequencies r2i,2 are different as well. 
Therefore, the Hamiltonian of the system of two NV spins 
interacting with the deformed cantilever reads 

H =Ho + V, 

Ho =^]4iwi(jI + ifiij2cr2 -I- fujjr [o^ a + 1/2) , (8) 

V =Ai (a+ + a)]^h [cos Oial + sin 6*1 (cr^ + o-f)] + 

A2 (a+ + a) ift [cos02O'2 + sin 6*2 (cr^ + ''■2^)] ■ 

Here the term huj-^ia^a -f 1/2) describes quantized os- 
cillations of the cantilever, tan0i_2 = and wi,2 = 
\l^\,2 + thereby denotes rii.2 the Rabi frequency of 
the transition between up and down spin states and 5 is 
the detuning between the microwave frequency and the 
intrinsic frequency of the spins. In the symmetric case 
hi = h2 = h, uji ^ UJ2 ~ ^ and Ai = A2 = A the Hamilto- 
nian (HI) recovers the previously studied modclSi^. The 
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coupling constants between the spins and the cantilever 
has the form Ai_2 = 5sA''bG'™20o, where gs ~ 2, /j,b is the 
Bohr magneton and G5"2 = ■jl-^tipl is the magnetic field 
gradient produced by the magnetic tips. Note that if the 
asymmetry of the deformation of the cantilever is strong 
Aft,(i?o) ^ 1, the distance between the magnetic tip and 
the nearest to the cantilever adjacent NV spin becomes 
very small z = hi leading to the very large magnetic 
field gradient G"° = 7^|-Btip| and to the large interac- 
tion constant Ai « gs^J,BG'^20■o/hl, while the coupling to 
the second spin becomes weak A2 ~ ffs/^BG™2ao/'i2 and 
therefore the relation ^ = ^ ^ 1 holds. This means 
that one can easily control the interaction between the 
magnetic tips and the NV spins simply by tuning the am- 
plitude of the external magnetic field and thus controlling 
the deflection of the cantilever AhiBo). In the symmetric 
case z = /ii,2 = h,Ah = 0, G5" = G^" = G^ ~ 10^ [T/m] 
and realistic values of the parameters are: ft, w 25 nm, 
A/(27r) w 0.1 [MHz], u},/{2n) w 5 MHz. In what follows 
the interaction constant A and the detuning between the 
cantilever frequency and the spin splitting A = lUi- — uj, 
A « 2A defines the time scale of the problem. A de- 
viation of the values of the constants from the values 
corresponding to the symmetric case reads 

Ai,2= A (fti,2) = A ± AA, Ai,2 = A (^1,2) = A ± Aw, 

I q2 

wi,2= \J{^± Anf + (52 w ± Aw, Aw = — {Ah) (9) 

CO 

AA= XAh, An = nAh. 

Estimates of the asymmetry parameter A/i(_Bo) describ- 
ing a deformation of the cantilever, for realistic materials 
and magnetic fields, follow in the next section. 



III. MAGNETIC-FIELD-INDUCED 
CANTILEVER DEFLECTION 

We propose to deposit a magnetic film on the can- 
tilever to explore alternatively magnetoelastic stress or 
magnetic-field induced torque to induce a controlled can- 
tilever deflection due to an external magnetic field Bq. In 
the first case we exploit the tendency of a film to develop 
magnetoelastic stress upon magnetization, and this stress 
induces a curvature of the thin cantilever substrate^"—. 
In the second case we exploit the magnetic torque^ 
T = rhxBo, where m is the total magnetic moment of the 
film, which is prepared to be oriented along the cantilever 
length (x-direction), and Bq is the applied magnetic field, 
oriented perpendicularly to the cantilever long axis (z- 
direction). In this geometry a cantilever deflection along 
the z-direction results. 

For a quantitative determination of the resulting can- 
tilever deflection we specify the cantilever dimensions as 
follows: length L = 3000 nm, width w — 300 nm, thick- 
ness tg ~ 30 nm. These parameters are well suited for the 
nano-fabrication of Si cantilevers^. They also represent 



a valid scenario of stress-induced free two-dimensional 
bending due to the large length-to-widtli ratio^. 

We assume that the cantilever is fabricated out of 
Si(lOO). We use the corresponding^ Young modulus Y = 
130 GPa and Poisson ratio 1/ = 0.279, and density p = 
2.33 X 10'^ kg/m~^. Such a cantilever has a mass of m = 
6.29 X 10~^^ kg. Its first three resonance frequencies^ 
are calculated from /^s = tsl3^{Y/{3p) f-^ /{4nL^) with 
13 = (1.8751,4.6941,7.8548) as 4.02, 25.2 and 70.6 MHz. 
This cantilever has a negligible deflection at its end due 
to its own weight of3pL^/{2Ytl) = 2.4 x 10"^^ m, which 
is more than six orders of magnitude smaller than the 
field induced defiection, as described next. 



1. Magnetoelastic-stress-induced cantilever deflection 

Magnetoelastic stress is responsible for the change of 
length of a bulk sample upon magnetization, and the re- 
sulting strain is known as magnetostriction^^i^. In films, 
a change of length of the film upon magnetization is not 
possible due to the bonding to the substrate, and the film 
develops a magnetoelastic stress. This stress induces a 
curvature of the substrate, which we exploit to dcfiect the 
end of the cantilever. The role of the external magnetic 
field is to induce a reorientation of the magnetization of 
the film from an in-plane (external magnetic field off) 
to an out-of-plane direction (external magnetic field on 
along z-direction). Such a reorientation of the magne- 
tization direction of the film with thickness tf induces a 
corresponding change of magnetoelastic stress. It induces 
a defiection of the free cantilever end given byi^ 

defUc = (10) 

For a magnetization reversal along the axes of a cu- 
bic system the magnetoelastic coupling coefficient Bi 
enters^^. For definiteness, we assume a Ni film of thick- 
ness 10 nm and we take Bi of bulk Ni, 9.38 MJ/m^. 
We note that the effective magnetoelastic coupling in 
thin films may deviate from its bulk value^^i^i^, but 
we stick to the bulk value for this proof of principle case 
study. With these assumptions we get a deflection of 
deflmc = 27.7 nm. A Ni film with bulk properties de- 
posited on the top surface of the cantilever has a tendency 
to contract along the magnetization direction. Thus, 
the cantilever would be curved upwards for zero external 
field (magnetization in-plane along x) , and it would curve 
downward for magnetic field on (magnetization along z). 
This is irrespective of the sign (along +z or along —z) 
of the external field. For AC magnetic fields (wmag) the 
deflection will change with 2 Wmag- 

The required defiection can be adjusted by varying 
the experimental parameters tf,ts,L, accordingly. Note 
that the magnitude of the external magnetic field does 
not enter. The only requirement is that it is large 
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enough to induce a magnetization reorientation. The ef- 
fective magnetic anisotropy of the magnetic film should 
be fairly small to achieve this. It can be tuned by adjust- 
ing film thickness, morphology, interface modifications, 
multilayer structures, and film composition to achieve 
this goal. The maximum field is given by the require- 
ment /iB-Bo < a-s pointed out above. This gives 
Bo < 32.7 mT. 



2. Magnetic-torque-induced cantilever deflection 

The deposition of a magnetic film on the cantilever 
opens also the possibility to exploit the torque induced 
by the external magnetic field along the z-direction on 
the magnetic moments oriented along the x-direction of 
the film to induce a cantilever deflection^. The magnetic 
torque T = rhx Bq is proportional to both the total mag- 
netic moment of the film m and the external magnetic 
field Bq. Here, we assume that the magnetic anisotropy 
of the film is large enough to ensure an in-plane mag- 
netic moment in presence of the out-of-plane field. Thus, 
the requirement on the magnetic anisotropy of the film 
differs from that discussed above for the magnetoelastic- 
stress-induccd defiection, as here an effective magnetic 
anisotropy favoring in-plane magnetization is required. 

For definiteness we assume a magnetic moment of 2fiB 
per film atom, and we refer to the atomic volume of 
bulk Ni, Patomic = 1.096 X 10~^^ m~^ to calculate the 
total number of Ni atoms in the film. We obtain for the 
torque- induced deflection^ rfe/Ztorquc = ATL? / {Ywt^), 
and from this we find 



deft 



torque 



/^atomic^^s 



(11) 



We obtain a defiection of 5.2 nm for the parameters 
quoted above at a film thickness of if = 10 nm for an 
external fieldi?o = 10 mT. Note that here the defiection is 
proportional to the external magnetic field, which allows 
a continuous control of the deflection. The magnetic film 
could be deposited on both sides of the cantilever, which 
would double the defiection in a given external magnetic 
field. Larger deflections at a given field are obtained by 
e.g. increasing the film thickness. 

We foresee that the exploitation of the magnetoelastic 
stress for achieving a well defined cantilever deflection 
is experimentally more challenging as compared to the 
exploitation of torque. One reason is the required exact 
tuning of the magnetic anisotropy. The torque approach 
is much more robust in that aspect, as only a sufficiently 
large energy barrier against out-of-plane magnetization 
is needed. Already the shape anisotropy of out-of-plane 
magnetization fulfills this requirement^. Therefore we 
focus now on the torque-induced deflection. 

To quantify the indirect interaction between the NV 
spins mediated by the interaction with the magnetic tip 
we consider again the asymmetry parameter A/i(i3o) = 



2\hi — h2\/{hi -\- /12). We calculate the magnetic tip-NV 
spin distances /ii, /12 as hi = ho—defl and = ho-\-defl, 
where ho describes the symmetric case, i.e. the cantilever 
end is at the center position between both NV spins, 
which is realized for zero magnetic field for the torque- 
induced defiection. This gives the asymmetry parameter 



Ah{Bo 



def Itoiquc _ IGtffiBBoL^ 



ho 



hoPntomicYtf 



(12) 



IV. EFFECTIVE HAMILTONIAN: INDIRECT 
INTERACTION BETWEEN SPINS 

Interaction of the NV spins with the magnetic tips 
leads to an indirect interaction between the NV spins. 
The Hamiltonian of the indirect interaction between the 
NV spins can be evaluated using the Frohlich method^ 



dt' [V (f) ,V m, V{t)^e- 



-iHoty^iHot 



(13) 

Taking into account eqs. (8) and in the rotating 
wave approximation we deduce 



cff 



1 




X2 




|sin2 




4 






1 






4 


sin Oi 


sin 62 


a 




Ai = 



Al A2 ^ A1A2 

ijJr — 1^1 , A2 = iO' 



'A2 



- UJ2- 



(14) 



l±\Ah 



Taking into account that sin0i_2 

1 and, therefore, sin^i 2 « — = sin6' we can rewrite the 
interaction constant in terms of the asymmetry parame- 
ter 



A1A2 sin 




1 - (Ahy 



(15) 



From cq. (|15|) we see that the asymmetry can enhance 
the interaction between the spins if 



a;2A2 



In the opposite case. 



> 1. 



< 1, 



(16) 



(17) 



the asymmetry lowers the interaction strength. Since 
the variance of parameters is relatively Iarge2ii2^ |^ « 
1 ^ 5 MHz, n w 0.1 ^ 10 MHz, 5 w 0.01 ^ 0.1 MHz, both 
cases given by eqs. ((T5)) and ((T7)) can be realized. 
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In particular from ea. p5|) we see that the interaction 
strength between spins depends on the asymmetry A/i 
and the dimensionless parameter a = Considering 
standard values of the parameters21i2^ w « fi, A ~ ^ 
from cq. (|15p for the interaction strength between spins 

we get: A[l + (a^ — l)(A/i)^]. We see that the depen- 
dence of the interaction strength between the NV spins 
on the asymmetry parameter Ah is not trivial. For small 
values of the parameter a < 1 the asymmetry arising 
from the deformation of cantilever leads to the reduction 
of interaction strength, while for a > 1 the asymme- 
try increases the interaction strength. In particular, for 
the values a = 3, Aft, = 0.5 the interaction between NV 
spins is three times larger than the interaction strength in 
the symmetric case ^ 3A. Stronger interaction between 
spins means larger entanglement. Therefore controlling 
the spin coupling strength we can influence the degree of 
entanglement as well. From the physical point of view 
small values of the parameter a sa ^ < 1 corresponds 
to a large detuning between oscillation frequency of the 
cantilever Ur and the Rabi frequency fls, while for a > 1 
we have the opposite case. 



INFLUENCE OF ASYMMETRY ON THE 
ENTANGLEMENT DEGREE 



In order to study the influence of the asymmetry on the 
degree of entanglement, we directly solve the Schrodinger 
equation corresponding to the Hamiltonian (|14p 



* dt 



(18) 



IV') = C'lit) |$+) + CaW 1$-) + Csit) |^+) + Ciit) . 



with the following notations 



2uj^ ' ujA 



B- 



1 (AhY 





1 - 




2w2A 


[jJa 


2uj ^ 


+ 1) 


IJa ' 



1 A/i, 



(20) 



Taking into account eqs. (18), (19) we can quantify the 
entanglement via the following expression 

ci\^m = \{r{t)\ay®ay\^pm- m) 

As a result, after straightforward but laborious calcula- 
tions for concurrence we deduce 

\C{t)\ = \ {C^{Q) ~ C^{Q)) cos2 At 

DC3{Q) + {F + B)Cm]^\ 



+e 



-2iFt 



+e 



-2iFt 



2iFt 



2iFt 



4F2 / 

[D^C3{Q)-D{F + B)C4{0)]- 
\ 4F2 {F + Bf 

' [DCM + {F^B)Ci{())f \ 
' [D^Cm + D [F - B)Cm]'^ 



(22) 



\ 4,F^{F-B) J 

l-DCm + {F + B) CijO)] [DCsjO) + (F - B) ^(0)] 
2F2 

[DCajO) -D(F + B) C^miD^CsiO) + D {F - B) ^(O)] 
2i^2 _ 5)2 



Here |<i>='=) and \ip^) are Bell states^"—. Taking into 
account eqs. dH]), (18) for the resolution coefficients we 
obtain 



Ci{t) =Ci (0) COS (At) - C2 (0) sin(At), 
C2{t) =C2 (0) cos {At) - Ci (0) sin(At), 



Ci{t) 



D'C3iO)-DiF + B)C^iO) 
2{F + B)F ^ 
D^Cs (0) + DiF-B)Ci (0) 
2{F-B)F 
-DCs{0) + {F + B)CiiO) 

= 2F ' 

, DCsiO) + {F-B)CiiO) 
^ 2F ' 



(19) 



VI. RESULTS 

Following eq. (22) the concurrence is plotted in Figs. 
13111] and [S] for a variation of two parameters a and Ah. 
The contour plot of \C{t)\ is presented in Fig. [S) 

From Fig. |4]we see that for small values of the param- 
eter a = 0, 25, the increase of the asymmetry Ah leads to 
a smaller concurrence. The explanation is that for a < 1 
the increase of the asymmetry Ah reduces the coupling 
strength between spins see eq. (12), while as we sec from 
Figs. |3l[5]and[6l for large a > 1 concurrence is increased 
with asymmetry. In particular, the contour plot of Fig. 
ini defines domains of maximal and minimal concurrence 
as a function of parameters Ah, a. For a < 0, 5 concur- 
rence is maximal for small asymmetry Ah < 0, 25, while 
for a < 0, 5 concurrence is maximal for large asymmetry 
Ah > 0,3. 
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reduces the interaction strength. In addition we found 



(1=0.25, Ah=0.4 



ty=0.50, Ah=0.4 



a=2.m, Ah=0.4 



15 20 25 30 



Time t 



FIG. 3. Time evolution of the concurrence |C(t)| for the 
given values of a and Ah. Other parameters read: n = 1, 
Ci = C2 = Cs = Ca — 1/2. Time-scale corresponds to the 
microsecond. 




a=l.5, Ah=0.35 



ff=1.5, Ah=0.45 



ff=1.5, Ah=0.55 



5 10 15 20 25 30 



Time t 



FIG. 5. Time evolution of the concurrence |C(t)| for the 
given values of a and Ah. Other parameters read: n — 1, 
Ci — C2 = Cs = Ci = 1/2. Time-scale corresponds to the 
microsecond. 




fy=0.25, Ah=0.1 



ff=0.25, Ah=0.3 



ff=0.25, Ah=0.5 



Time t 



10 15 20 25 30 



FIG. 4. Time evolution of the concurrence |C(t)| for the 
given values of a and Ah. Other parameters read: n = 1, 
Ci = C2 = Cs = Ci — 1/2. Time-scale corresponds to the 
microsecond. 



"■5 Ah, [d.u.] 

0.35 



a, [d.u.] 



1.5 



VII. CONCLUSIONS 



FIG. 6. Contour plot of the concurrence |C(t)| for t = 13.3. 
Other parameters are: n = 1, Ci = C2 = C3 = C4 = 1/2. 



One of the main challenges for NV spin-based nano- 
elcctromcchanical resonator has been the achievement of 
a high controlled degree of entanglement and a strong 
coupling between NV spins. With this in mind, we pro- 
posed in this work a new type of nano-electromechanical 
resonator, the functionality of which is based on the NV 
spin controlled by an external magnetic field. In particu- 
lar, we suggest to deposit a thin magnetic Ni film on the 
Si(lOO) cantilever to exploit alternatively magnctoelas- 
tic stress or magnetic-field induced torque for inducing a 
controlled cantilever deflection upon acting with an ex- 
ternal magnetic field. We have shown that, depending on 
the values of parameter a = ^{^2 ■, the induced asymme- 
try of the cantilever deflection substantially modifies the 
characteristics of the system. In particular we demon- 
strated that is if a > 1 the asymmetry enhances the 
strength of the interaction between the NV spins at least 
three times ^ 3A, where A = lOOkHz is the maximal 
coupling strength between NV spins for the symmetric 
model reported in2i. However, for a < 1 the asymmetry 



that for a > 1 entanglement is maximum for the case of 
a large asymmetry Aft,, while for a < 1 the entanglement 
is maximal in the small asymmetry case (See Fig. 6). 
The values of the parameter a can be changed efficiently 
via the change of the detuning between the oscillation fre- 
quency of the cantilever and the spin splitting frequency 
A = — oj. This can be used as an effective tool for a 
practical implementation of our theoretical proposal for 
controlling the entanglement and the interaction strength 
between NV spins in the experiment. 
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